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It is proposed that the newly discovered X(3872) is a J^*^ = l^""" D^'D''* hadronic resonance 
stabilized by admixtures of ujj /ijj and pJ jii). A specific model of the state is constructed and tests 
of its internal structure are suggested via the predicted decay modes D^D^tv", D'^D'^'J, n'^w~ J/tp, 
and tt'^tt^ TT° J/tli. 
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' The Belle collaboration has recently announcedQ the discovery of a resonance, X(3872), in the tt+tt J/ip subsystem 
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This state, which has been confirmed by the CDF collaboration has attracted some attention because of its 
unusual properties. Specifically, the state appears to be too heavy to be a ID charmonium state and too light to be 
2P charmonium or a cc hybrid. See Ref. 3;] for a detailed assessment of possible charmonium assignments and decay 
T-H modes. 

Alternatively, the proximity of the state to DD* threshold strongly suggests that the X may be a weakly bound 
I DD* resonance 0, 0,1^ 0, Hi , sometimes called a mesonic 'molecule' or a 'deuson' 9]. This is an old idea which has 
I been applied to a variety of mesons with unusual characteristics such as the V'(4040jjlflj|, /i (1420) [Hill, »C1440)[T^. 
Oh! /j(1720)[il,ao(980), and/o(975)[ll[l3. 

JL , In this note I assume that the AT (3872) is indeed a DD* resonance and present a detailed analysis of its expected 
properties based on a simple model of quark interactions. This model incorporates the nonrelativistic quark model 
(~| , with additional dynamics due to pion exchange. The idea is to capture the predictive power of a microscopic formalism 
I!" ' of short range quark dynamics along with the important long range dynamics mediated by pion exchange processes. 
Versions of this idea have been applied to baryon-baryon interactions since the 1980'sT6j, where, of course, pion 
exchange is of fundamental importance; another variant has recently enjoyed some vogue in baryon physics 17|. 
^1 It is natural to expect that the putative DD* bound state is in a relative S-wave since this is typically where 
inter-hadron forces are strongest. In this case pion-mediated interactions (see below) favour the isoscalar channel, 
which in turn implies a J^'~' = 1++ state. I will therefore henceforth refer to the bound state interpretation of the 
AT (3872) as the Xci(3872). The remainder of this paper focusses on the properties of this state. 

Although pion exchange forces dominate the structure of the Xci (in analogy to the deuteron), short range quark 
dynamics are present and assist in binding the Xci via mixing to hidden charm vector- J/^/; states. Indeed, ojJ/ip and 
pJ /ip are very nearly degenerate with DD* and one must expect some admixture of these states - an effect which 
will be strongly enhanced by the near-zero energy denominator. Such mixing is also important in driving possible 
decay modes of the Xci and is therefore central to determining its properties. Finally, the binding energy of the Xci 
is comparable to mass differences in the available charge channels and one can expect strong isospin violating effects 
in this resonance. This heretofore unexplored dynamics is thoroughly examined in the following. Detailed predictions 
of binding energies and branching fractions are presented along with possible experimental tests of Xci structure. 



II. DD* DYNAMICS 



Long range pion exchange effects are expected to dominate the physics of a weakly bound state such as the Xci- 
Nevertheless, as discussed above, short range quark interactions can give rise to important mixing effects. We therefore 
consider a model which appends pion exchange dynamics to the nonrelativistic quark model. The model is used to 
extract effective interactions for DD* — DD* , DD* — loJ /if), and DD* — loJ jip scattering. These interactions are then 
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employed in a nonrelativistic coupled channel Schrodinger equation to extract bound state properties (one expects 
the nonrelativistic formalism to be accurate for weakly bound states of relatively massive components as is the case 
with the Xci)- 



A. Quark Exchange Induced Effective Interaction 



The quark model employed here assumes nonrelativistic quark dynamics mediated by an instantaneous confining 
interaction and a short range spin-dependent interaction motivated by one gluon exchange. The colour structure is 
taken to be the quadratic form of perturbation theory. This is an important assumption for multiquark dynamics 
which has received support from recent lattice computations for both confinement 18] and multiquark interactions 
The final form of the interaction is thus taken to be 
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where A is a colour Gell-Mann matrix, as is the strong coupling constant, h is the string tension, rrii and rrij are the 
interacting quark or antiquark masses, and cr is a range parameter in a regulated spin-spin hyperfine interaction. The 
parameters used were as = 0.59, h = 0.162 GeV^, a — 0.9 GeV, and 0.335, 0.55, and 1.6 GeV for up, strange, and 
charm quark masses respectively. Relevant meson masses obtained from this model are p = 0.773 GeV, J /ij) = 3.076 
GeV, D = 1.869 GeV, and D* — 2.018 GeV, in good agreement with experiment. 

Meson-meson interactions are obtained by computing the Born order scattering amplitude for a given process [13. l20j . 
Because of the colour factors in Eq. lathis amplitude necessarily involves an exchange of quarks between the interacting 
mesons. Thus the leading order DD* interaction couples DD* with hidden charm states such as pJ/'0 and ujJhb. 
This amplitude may be unitarized by extracting an effective potential and iterating it in a Schrodinger eauation|l2|. 
The method has been successfully applied to a variety of processes such as KN scattering '2ll| and J/-0 reactions 
relevant to RHIC physicSj22j. It has even proven surprisingly useful for relativistic (and chiral) reactions such as tttt 
scattering(]J,|20|. 

The S-wave Born order scattering amplitude for DD* — loJ/iI) scattering is shown in Fig. ^ Here DD* refers to the 
isoscalar positive charge parity state l/^/2{DD* + DD*)^. The scattering amplitude is dominated by the confinement 
interaction of Eq. (this is in contrast to light meson scattering which is dominated by the hyperfine interaction) . An 
effective potential is extracted by equating the scattering amplitude to that obtained for point-like mesons interacting 
via an arbitrary S-wave potential. It is convenient to parameterize this potential as a sum of gaussians: 



n^^a.e-^/2^1 (4) 

i 

The fit to the quark level amplitude is illustrated in Fig. ^ (left panel) and the resulting potential is shown in Fig. 
n (right panel). The distinctive "mermaid potential" seen here is due to destructive interference between diagrams in 
the quark level amplitude. Thus details of the potential are sensitive to the assumed microscopic interaction, however, 
its general shape and strength are quite robust [l^ . 

The derived parameters of Eq. Hwere ai = 6.35 GeV, h = 1.166 GeV'^ and 02 = -6.82 GeV, &2 = 1-096 GeV'^. 
These parameters were obtained for the confinement portion of the "prior" form of the scattering amplitude. The 
"post" form yields ai = 3.82 GeV, h = 1.20 GeV^^, 02 = -4.21 GcV, and 62 = 1-125 GeV^^. Post and prior forms 
of a scattering amplitude refer to different schemes for constructing the time evolution operator which exist in the 
scattering of composite systems. In principle these give rise to the same scattering amplitude, but approximations 
and inaccurate wavefunctions can cause slight differences as indicated in the figure. We employ the average potential 
indicated by the solid line of Fig. ^ (right) in the following. Finally, the isovector DD*- pJ/ip effective potential is 
identical to its isoscalar analogue. 

The mermaid form of the effective potential implies that quark exchange effects can cause binding in the coupled 
DD*, ujj/ip or pJ/ip systems; however, direct computations indicate that the potential depth is not sufficient to form 
a resonance. We therefore turn to an examination of pion exchange induced dynamics in the DD* system. 



B. Pion Exchange Induced Effective Interaction 



I choose to follow the method of Tornqvist [1^] in constructing an effective pion- induced interaction. This is based 
on a microscopic quark-pion interaction familiar from nuclear physics: 
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Z = J d^x4,{x)ri5T''i'{x)d^^'^{x). (5) 

Here — 92 MeV is the pion decay constant, r is an SU(2) flavour generator, and g is a coupling to be determined. 
The eflFective potential is derived by projecting the quark level interactions onto hadronic states in the nonrelativistic 
limit. In the case of pseudoscalar- vector states one obtains [T^ 



where 



(6) 



and 



C{r) = ^^, (7) 

rM.C(,.)(i + l: + ^), (8) 

yo^^|-1.3MeV. (9) 

The matrix elements refer to S- and D-wave components of the pseudoscalar- vector state in analogy with the deuteron. 
The strength of the interaction has been fixed by comparing to the ttNN coupling constant via the relationship 
5^AfAr/47r = 25/18 • m'^g'^/f^. This allows a prediction of the D* decay width which is in good agreement with 
exDeriment|l3| . The parameter /i is typically the pion mass, however, one can incorporate recoil effects in the 
potential by setting fi^ — — [mv — mps)"^ ■ The results presented here are insensitive to the value of ^ and I take 
/J, = 130 MeV in the following. Finally, the coupling 7 is a spin-flavour matrix element which takes on the following 
values: 7 = 3 for / = 0, C = -)-; 7 = 1 for J = 1, C = -; 7 = -1 for / = 1, C = -h; and 7 = -3 for / = 0, C = -. 
Thus the isoscalar positive charge parity channel is the most likely to form bound states and subsequent discussion 
focusses on it. 

The potential of Eq. [S] is an illegal quantum mechanical operator and must be regulated, typically with a dipole 
form factor. The regulator scale, A may be fixed by comparison with nuclear physics; for example NN interactions 
yield preferred values for A in the range 0.8 GeV to 1.5 GeV depending on model details. Alternatively, reproducing 
the deuteron binding energy requires A « 0.8 GeV. Tornqvist has employed an intermediate value of A = 1.2 GeV 
which is appropriate for D mesons and this is taken as the canonical cutoff in the following. 

Integrating the coupled S/D wave system for the 1++ BB* system yields a bound state of mass 10562 MeV, 
in agreement with Ref. [13. Similarly a + BB* bound state of mass 10545 MeV arises from this formalism. 
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Unfortunately, D mesons are sufficiently light that the DD* system does not bind with canonical parameters. However, 
the combined pion and quark induced effective interactions are sufficient to cause binding. The properties of this bound 
state are explored in the next section. 

III. PROPERTIES OF THE Xci(3872) 

The proceeding considerations indicate that the isoscalar positive charge conjugation sector is the most likely to 
bind in the DD* system. Furthermore the small branching fraction of the Xci to TnrJ/ip implies a small isovector pJ/ip 
component in the Xci wavefunction. Thus a good initial study is provided by the coupled channel l/^/2{DD* +DD*)g, 
l/\/2{DD* + DD*)^jj, and cuJ/ip system. Utilizing the potentials of Eqs. Qla-ndEland meson masses of D — 1.869 
GeV, D* = 2.01 GeV, w = 0.78 GeV, and J/V' = 3.1 GeV yields a single bound state of mass 3.872 GeV without 
adjusting any parameters, in remarkable agreement with the mass of the X. The Xci wavefunction is plotted in Fig. 
121 (left panel); one sees typical deuteron-like wavefunctions with strong D-wave and ajj/i/; components. All three 
components are required to achieve binding for this state. 

Although encouraging, this result must not be taken too seriously because the binding energy is comparable to 
the difference in energies of the various relevant charge channels. Thus isospin breaking effects are expected to be 
important and must be incorporated in the formalism. This is achieved by including isovector channels and allowing 
for differing thresholds. Restricting attention to nearby vector meson - J/tp states and neglecting the coupling to 
charmonium states^J yields six possible channels: 

• pJ/tP at 3.8679 GeV 

• -1={D°D°* + D"D°*)s at 3.8712 GeV 

• -^{D°D°* + D^D^*)d at 3.8712 GeV 

• ^(L>+i:)-* + i:'-L>+*)s at 3.8793 GeV 

• ^{D+D-* + D-D+*)d at 3.8793 GeV 

• LoJ/iP at 3.8795 GeV. 

One sees an immediate problem: the threshold for pJ /ip is too low to allow a resonance at 3.872 ± 1 GeV. However, 
the mass of the p is rather poorly defined due to its large width and some leeway in fixing threshold for this channel 
is permissible. I therefore adopt the simple prescription of setting the p mass equal to that of the u at 0.7826 GeV. 
Varying this prescription made negligible changes to the following results. The quark level coupling of {DD*)d states 
to S-wave p- or lo-J /ip states is small and is neglected. The final step is to form effective interactions from appropriate 
combinations of isospin basis interactions. For example the D^D^* — D^D^* interaction is given by Eq. Elwith 7 = 1. 
The resulting numerical six channel problem must be studied with some care because the binding energies are small 
relative to the natural scales of the system. 

It is of interest to study the properties of possible bound states as a function of their binding energy. This has been 
achieved by allowing the regulator scale to vary between 1.2 and 2.3 GeV. Binding is seen to occur for A larger than 
approximately 1.45 GeV. Wavefunction coefficients (defined as / Iv'aP where a is a channel index) are shown as a 
function of binding energy in Fig. [5] (right panel). It is clear that the D^D^* component dominates the wavefunction, 
especially near threshold. However, the D~^D~* component rises rapidly in strength with isospin symmetry being 
recovered at surprisingly small binding energies (on the order of 30 MeV). Alternatively, the ujj/ip component peaks 
at roughly 17% at Eb ~ 9 MeV. The contribution of the pJ /ip wavefunction remains small, peaking at less than 1% 
very close to threshold. 

It is possible to estimate decay rates in a simple fashion once an explicit wavefunction is known. This is because 
the small binding energy of the Xci implies that its constituent particles are nearly on-shell and therefore Xci decay 
amplitudes are well approximated by constituent decay amplitudes. Thus, for example, the tt+tt" J/i/; decay mode 
arises predominantly from the pJ/ip wavefunction component. The channel strengths of Fig. |2] therefore allow simple 
estimates of a variety of branching fractions based on the widths of the D*, lo, and p (decays of the D and Jf-p 
mesons are neglected here but can be computed easily). The results for a variety of modes are presented in Table I 
as a function of the binding energy. 

The broadest particle in the Xci system is the p with a width of 150 MeV and it is the pJ /ip component which 
has the largest branching fraction, even though it is strongly suppressed in the wavefunction. The next strongest 
mode is provided by the tt^tt^tt^ decay of the uj which is enhanced relative to other modes due to strong mixing 
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r(1/GeV) EB(MeV) 



FIG. 2: (left) Three Channel Isoscalar Wavefunction Components, (right) Component Strength vs. Binding Energy. 

with LoJ /ij:. Unfortunately, only a very rough upper limit on the total width of the D^* exists|23| so estimates of the 
D^D^TT^ and D'^D'^j decay widths are essentially useless. The figures in the table have been obtained by assuming 
that r{D°* D°j) « 25 keV and r(L»°* ^ D°tt°) « 43 keV; both of these estimates are anchored in D^* decays 
and should be reliable. Notice that the D^n^ mode is closed. All other possible decay modes of the Xci are relatively 
small, although the ■k^^J/iI) mode may be of interest if it is detectable. 









TABLE I: 


Some Decay Modes of the Xci(3872) (keV). 
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IV. CONCLUSIONS 

I have argued that the A(3872) is a J^'~^ = 1++ DD* hadronic resonance with important admixtures of pJ/ip 
and LoJ/il^ states, dubbed the Xci- This assertion is supported by detailed computations in a microscopic model 
which incorporates pion and quark exchange interactions. The model has been heavily tested on nuclear physics and 
meson-meson scattering data and can be regarded as reasonably reliable. The 1++ Xci is the only DD* state which 
binds; no other J^*^ or charge modes exist in this model. Furthermore, no DD molecules are expected. It is likely, 
however, that a rich D*D*, BB* and B*B* spectrum exists. Thus the discovery of the A(3872) may be the entree 
into a new regime of hadronic physics which will offer important insight into the workings of strong QCD and should 
help clarify many open issues in light quark spectroscopy. Indeed, the experimental and theoretical analysis of heavy 
molecules is simplified because of their hidden flavour components. 

It is clear that further experimental studies of the A (3872) are of great importance. For example, determining its 
spin and parity are of immediate concern. The fact that the A is polarized in i? — )■ KX will help greatly in this. 
Furthermore, detecting a tt'^tt^J/tP decay mode would immediately eliminate the Xci interpretation of the A. 

It is also important to gather enough events to reconstruct the invariant mass of various subsystems such as tt+tt^ in 
7r+7r~ J/f/; (which should peak at the p mass). Perhaps a more interesting test would be the invariant mass distribution 
of the Tr+TT^Tr" subsystem in the 7r+7r~7r° J/V' decay mode, which should have all of its events near the edge of phase 
space due to the narrow width of the virtual lo. It is therefore encouraging that the SirJ/tp decay mode is roughly 
1/2 the strength of the 2nJ /ip mode. Although some events will be lost due to the decreased efficiency in detecting 
neutral pions, this deficit should be made up by the new data being collected at Belle and BaBar. 

Lastly, determining branching fractions, especially those arising from different wavefunction components such as 
D^D^TT^, TT^TT^ J/ip, and Tr+vr^vr*' J/i/), would help greatly in pinning down the internal structure of the A and provide 
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an intriguing glimpse into a new realm of hadronic physics. 
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Coupling of the Xci to charmonium states is negligible because the wavefunction squared at the origin scales as \/ finD" Eb 
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by S-wave quark pairs and possible charmonium states are all orbital excitations. 



